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Abstract. We use exact diagonalization and the modified Lanczos method to 
study the finite energy and finite momentum spectral weight of the longitudinal and 
transverse spin excitations of the anisotropic zig-zag ladder. We find that the spin 
excitations form continua of gapless or gapped spinous in the different regions of the 
phase diagram. The results obtained are consistent with a picture previously proposed 
that in the anisotropic case there is a transition from a gapped regime to a gapless 
regime, for small interchain coupling. In this regime we find a sharp low-energy peak 
in the structure function for the transverse spin excitations, consistent with a finite 
stiffness. 



Spin dynamics of the 5=1/2 antiferromagnetic zig-zag ladder with anisotropy 



2 



1. Introduction 

Recently |I|] it has been suggested that a two dimensional spin system, CS2CUCI4, 
has an excitation spectrum that can be described, similarly to the one dimensional 
case [0], by a continuum originated from pairs of spin 1/2 spinous. The standard one 
dimensional Heisenberg model is known to have fractional states where the usual spin 
1 magnons are replaced by pairs of deconfined spin 1/2 topological excitations called 
spinous The characteristic low excitation energy coherent peaks that appear in 
the dynamical susceptibility are in this case replaced by a continuum. This property 
has been verified experimentally for several quasi-one dimensional spin 1/2 systems 
like CPC 0], KCuF^ |^ and copper benzoate where a description in terms of a 
nearest-neighbor Heisenberg model is assumed to apply. 

Real materials are however neither strictly one-dimensional nor the interactions 
are of the simple Heisenberg nearest-neighbor form. Originally CS2CUCI4 was taken 
as a quasi-one- dimensional system ^ but a more careful estimate of the inter-chain 
parameters revealed that they are of the same order as the intrachain interactions |l[]. 
The interlayer coupling is estimated to be two orders of magnitude smaller implying 
that the system is essentially two-dimensional. It forms a triangular lattice in the 
copper {Cu^~^, S = 1/2) planes and constitutes therefore a frustrated system. 

Frustrated systems have attracted considerable interest. Using a large- bosonic 
expansion it has been predicted M that the presence of frustration may counteract 



the staggered fields responsible for confinement |10] and lead to deconfinement of the 



spinous. In a non-frustrated system it is known that the low energy modes are spin- 



1 magnons. The presence of frustration leads to a non-coUinear order parameter 
which can be parametrized by three real scalar numbers. The representation of the spin 
operators in terms of bosonic operators contains an hidden internal gauge symmetry 
under which the scalar link fields transform either as charge +2 or charge -|-1 scalars 
0. It has been predicted a long time ago that if the charge +2 scalars condense into a 
Higgs phase [|12| then the unit charges are not confined and the spinous remain free fTT 



However the Higgs phase is just one of the possible phases predicted to occur in the 
frustrated two-dimensional lattice. The CS2CUCI4 system provides a first experimental 
example of a two-dimensional system with a spectrum consistent with the existence of 
spinous. 

One may think of a simpler system like the zig-zag ladder, where only two chains 
are coupled, and study the excitation spectrum of this system as a first step towards 
understanding the two-dimensional triangular lattice obtained in the limit when several 
ladders are coupled to each other. In usual spin ladders an infinitesimal coupling between 
chains leads to a behavior qualitatively different from the Id case. Therefore it is 
interesting to see if in the case of the zig-zag ladder there is a qualitative difference from 
the single-chain case. It has been found before that if the next-nearest-neighbor (nun) 
interaction is small enough the system still behaves qualitatively as the single-chain case. 
Therefore it seems reasonable that for small couplings we might find similar features 
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characteristic of the single chain, in particular a spinon originated spectrum. The main 
question to be answered is what happens for large couplings where the role of the nnn 
coupling is important. We will show that indeed for all couplings a spinon description 
holds. A direct comparison with the experimental results for CS2CUCI4 would require 
a full 2d calculation but our results for the zig-zag ladder give a first indication that the 
nnn coupling does not lead to a coherent energy spectrum. 

In the context of a frustrated two-leg ladder similar to the zig-zag ladder 



it has been shown that the spinous survive as the elementary excitations in a 
spontaneously dimerized ground state but become massive. A local Z2 symmetry related 
to independent translations by one lattice spacing on each chain is spontaneously broken 
and leads to a non-vanishing dimerization for strong enough frustration. This symmetry 
breaking leads to kinks as elementary excitations which are massive. These kinks have 
been shown to have S = 1/2 and therefore at least two must be created. 

In general if a term removes explicitly the Z2 degeneracy betweeen the two 
dimerized ground states this leads to soliton confinement. An example is to introduce 
explicit dimerization in the Hamiltonian. The role of explicit dimerization has also 



been addressed in the context of spin-Peierls systems like CuGeOs [14|, NaV20, 



||T5| or C u{N 0-i) 22. hD20 ||T6[. The excitation spectrum of these systems is however 
considerably different from the chains without dimerization. The spectrum is gapped 
and the lowest energy excitations are coherent spin-1 magnon peaks [0 that separate 
from the continuum that appears at higher energies. The spectrum of these systems is 
actually closer to integer spin chains |T8| or to spin-ladders ||T^. The effects of explicit 
dimerization have received renewed attention recently . 

Another source of interest in the zig-zag ladder is that it has been proposed that 
in the anisotropic case incommensurate quasi-long-range spin correlations should be 
observed. Also a gapless chiral phase has been predicted to occur |21]. In this work 
we focus our attention on the combined effects of a next-nearest neighbor frustrating 
interaction and of anisotropy. In a two-dimensional non-frustrating system the spectrum 
is coherent and composed of magnons. The addition of frustrating terms may lead to 
deconfined spinous. In the zig-zag ladder when the nnn frustrating coupling is absent 
the system has a spectrum determined by gapless spinous (it is the case of the simple 
Heisenberg chain). Adding frustration it is expected that the spinous will remain 
deconfined. In the isotropic case where the effect of the nnn coupling is to dimerize 
the spinous are massive at sufficiently strong nnn coupling. In the anisotropic case 
at strong enough nnn couplings we expect the system to have a transition from the 
intermediate dimer phase to a gapless phase where we expect the deconfined spinous to 
be massless. 
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2. Hamiltonian 



The anisotropic zigzag ladder is defined by the Hamiltonian 



z 

+1 



The spin operators refer to spin 5=1/2 states, while the summation ^ = 1, runs 
along the "rib" of the zig-zag ladder. We shall parameterize the interactions by the 
coupling parameter j = J^^/J/*-^ and by the anisotropy parameter Jf/J^"^ = A = 
J|/J^^. (The isotropic case reduces to j = J2/J1 and A = 1.) We will set Jf^^ = 1 
henceforth. The nearest-neighbor Heisenberg chain with anisotropy corresponds to both 
the weak-coupling (Ji = 0) and to the strong-coupling (J2 = 0) limits of the zig-zag 
ladder. The spectrum is gapless for the case of XY anisotropy, |A| < 1, as shown by 
the Bethe ansatz p2[. The excitation spectrum consists of spin-1/2 particles dubbed 
spinous. Since flipping one spin represents a spin-1 excitation, the spinous can only be 
created in pairs. Therefore the conventional spin 1 magnons are deconfined into spin-1/2 
spinous that propagate incoherently. 

The isotropic case has been studied before E!], EBI, HTl, as a function 



of the coupling parameter, j = J2/ Ji- As j increases, the system goes from gapless 
(single chain) to a dimer phase and then to a spiral phase, where the structure factor 
has a maximum at a momentum n/2 < q < 7c. The system has a spin gap in these last 
two phases, and it therefore only displays short-range order. In the limit that the intra- 
chain interaction is much larger than the inter-chain interaction {j — >■ 00) the two chains 
decouple and a gapless single chain behavior is recovered. It has been argued that this 
only happens, strictly speaking, at j = 00: the spin gap becomes exponentially small 
as j grows, but it remains non- vanishing . Recently, on the other hand, it has been 
proposed that incommensurate quasi-long-range spin correlations should be observed if 



easy-plane (XY) anisotropy is included in the zigzag ladder This is argued to be 
due to the presence of a "twist" term that results from the inter-chain interaction. It has 
been proposed that there is one gapless mode and one mode with a gap in the regime 
of strong XY anisotropy in the inter-chain coupling. Another prediction of this work is 
the existence of spontaneous local spin currents. This, however, has been refuted in ref. 



29| . Also, other recent numerical work has failed to confirm the gapless nature of 



the groundstate in the anisotropic XY case at weak interchain coupling. Recent Density 
Matrix Renormalization Group (DMRG) results ||^ suggest, however, that the zig-zag 



ladder does indeed show a gapless chiral phase as predicted in ref. ^1 . 

Also, recently an analysis of the exact properties of such finite systems was carried 
out, looking at various correlation functions and the structure of the spectrum both in 
the isotropic and the anisotropic cases |3^. The spin stiffness of the zig-zag ladder was 
calculated, and it was found evidence for a gapless regime at weak coupling that survives 
the thermodynamic limit in the case of XY anisotropy. This was also concluded looking 
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at the level crossings to detect the phase transition between the two regimes p2[ using 



a previously proposed procedure to detect the dimer transition at strong interchain 



coupling |23|. The same method was also recently used in reference ||3l 

In this work we will study the structure function and the spectral weight of the 
spin excitations both for the longitudinal and the transverse correlations. Our results 
are consistent with previous conclusions that there is a transition to a gapless regime 
at weak coupling if anisotropy is present. The results indicate a continuum of gapless 
excitations in the transverse correlations in the XY case, and a continuum of gapped 
excitations in the isotropic case. 

3. Spectral weight 

The structure function is defined by the overlap of two states coupled either by the 
longitudinal or the transverse spin operator P], 

S,M. ^) = jrjll e^'"" / dte^^' < ^r(i)^r+i?(0) > (2) 

l,R •'-•^ 

where and v are cartesian components. At zero temperature we obtain therefore 

S^,{q, uj) = Y. M^S{u + Eg- Ex)6{q + h - h) (3) 

A 

where Eq is the groundstate energy, Ex is the energy of an excited state, u is 
the excitation energy (energy difference to the groundstate) and q is its momentum 
(momentum difference to the momentum of the groundstate /cq) and the spectral weight 
is defined by 

= 2n\< G\S^{q)\X>f . (4) 

where S'^{q) is the Fourier transform of the spin operator. We will calculate the structure 
functions S+-{q,uj) and Szz{q,^) which probe the transverse and the longitudinal spin 
excitations, respectively. 

The single chain case was studied before both by Bethe ansatz and using 
numerical diagonalization of small systems 0. The continuum of excitations is 
contained in the thermodynamic limit between two lines: the bottom one is the single- 
spinon dispersion, ei{q), and the upper one is the maximum energy resulting from the 
combined effect of two spinous, e^iq). In the isotropic case (A = 1) the lines are defined 
by 



Jl 2 

en(g) 



TT 



Jl 

and in the XY case (A = 0) are defined by 



sin(g)| 

(5) 



sml- 



en(g) 



sin(g)| 

(6) 



sm,- 
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In the XY case the longitudinal structure function can be calculated exactly 
since this system is equivalent to free spinless fermions and it is given by ||^ 

^ e(a;-sin(g))9(2sin(g/2)-a;) 

Szz{q, ^) = 2 ; ^ = (7) 

/4sin2(g/2) -cj2 



In the isotropic case there is no exact solution but Miiller et al. ^ proposed an ansatz 
that fits very well both numerical results for small systems and various experimental 
results where a description in terms of a single chain is expected to hold. The Miiller 
ansatz is 

SUq, ^) = , ^ e{oo - ei{q))e{eu{q) - u) (8) 
^uj^-ei{q) 

where 9 is a step function and A a constant || . This function diverges at the lower 
boundary while in the XY case it diverges at the upper boundary [Q. At momentum 
TT the divergence is stronger and it diverges as Szz ~ 

In the thermodynamic limit the structure function eq. (3) can be written as a 
product 

S^^{q,u) = M>'{q,uj)D{q,u) (9) 

where M^{q, uj) is the continuum limit of the spectral weight originating in the overlaps 
eq. (4) and D{q,uj) is the density of states. In the isotropic and in the XY case the 
density of states is finite and nearly constant close to the low-energy threshold and it 
diverges at the upper threshold. On the other hand M'^{q, u) is constant in the XF-case 
and it diverges at the lower threshold in the Heisenberg case. The structure function as 
a consequence diverges in the lower threshold for the Heisenberg chain and it diverges 
in the upper threshold in the XF-case . For any finite system the density of states is 
a set of delta functions at the excitation energies. 

Using field theory it is also possible to determine the dependence of the transverse 
structure function close to the lower threshold. In the single chain case the transverse 
function is given by 

S+-iq,-^) - . , \. ^^y, ,J = 0, A = (10) 

and therefore 

However the finite energy structure function is not known analytically. 

The ladder case is more involved. We will use exact diagonalization of finite systems 
together with the modified Lanczos method IR7|. 



4. Numerical results 



Let us begin by recalling the quantum numbers of the groundstate as a function of 
the size N for the S* = 1/2 zig-zag antiferromagnet. Periodic boundary conditions 
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are imposed throughout. The groundstate is a spin singlet in general due to the 
antiferromagnetic interactions. The system has three well defined regimes: (a) strong- 
coupling, (b) intermediate coupling and (c) weak-coupling. Consider the isotropic case 
first. For strong enough coupling between chains, j = J2/J1 < 1/2, it has either 
momentum vr for N = An + 2 or momentum for N = An. For intermediate couplings 
(j > 1/2), on the other hand, the momentum oscillates between and tt as a function of 
the coupling parameter j and of the system size |^9[|. There are several points along 
j in this regime where the corresponding energy levels for these two momentum values 
cross. The groundstate is degenerate at these points, and this is refiected by peaks 
in the dimer correlation function ||3^, Such level crossings grow in number as the 
system size grows, and this indicates that the two singlet states in question are in fact 
degenerate in the thermodynamic limit. By the Lieb-Schultz-Mattis theorem, this is 
consistent with a spin gap in the excitation spectrum that survives the thermodynamic 
limit in the weak-coupling regime j > 1/2. 

The spectrum of the anisotropic S = 1/2 XX Z zig-zag ladder has also been 
studied previously in the strong-coupling regime up to the Majumdar- Ghosh line 
(0 < j < 1/2). ^2[ A gapless regime occurs for XY anisotropy A < 1 and strong 
coupling j < jci(A); an Ising antiferromagnet along the rib of the zig-zag that shows a 
spin gap in the excitation spectrum occurs for A > 1 and j < jci(A), and a dimer phase 
regime that also has a spin gap exists at j > jci(A) and any A. The line j = jci(A) 
separates the gapless phase from the dimer phase for A < 1, while it separates the dimer 
phase from the (Ising) Neel phase for A > 1. The line at A = 1 and j < jd separates 
the XY gapless phase from the Ising phase. 

It was found I^Sj that there is a transition from the gapped intermediate coupling 
regime to a weak-coupling gapless regime. In the intermediate coupling regime the 
two lowest states are two states with 5*^ = 0, of momenta = 0, tt. There is a line, 
jc2(A), where the first excited state becomes a. Sz = ±1, k = tt/2 state signalling the 
doubling of the periodicity and leading to a gapless regime as confirmed from the spin 
stiffness tensor highest eigenvalue [^. The curve jc2(A) shown in Fig. 4 of ref. 
separates a spin-gap (dimer) phase from a gapless phase at small interchain couplings. 
As expected, the value of jc2 grows near the isotropic point. (It should tend to j = cx) 
at A = 1 according to White and Affleck, |2^ but finite-size effects gave a finite value). 

Let us now analise the spectral weight eq. (4) at various points in the phase diagram 
parametrized by j and A. We will focus our attention on two classes of parameters. We 
will consider the isotropic case (A = 1) and the Xy-case (A = 0) varying in both cases 
the interchain coupling, j. 

Let us begin with the single chain case for both values of A and let us consider 
the particular case A^ = 16. In Fig. 1 we show the lowest energy levels (taking the 
groundstate as the zero of energy) for Sz = and 5*^ = 1 and a given momentum for the 
Heisenberg chain and the XY chain. In the Heisenberg case the states are organized 
into spin multiplets due to the SU (2) spin invariance. The groundstate is a spin singlet 
with momentum zero. The first excited state is a spin triplet with momentum k = n 
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and the next state is another spin singlet but with momentum k = n. In the A = 
case the first excited state is now a state with 5*^ = 1 and momentum k = tc and the 
next state is fourfold-degenerate with Sz = and momentum k = Air/N or momentum 



k = 71, OT Sy = 2 and the same momentum values 132 



The structure function Szz only couples the groundstate to states with 5*^ = and 
5*^ = 1- On the other hand only couples the groundstate to states with Sz = I 
and St = I {St is the total spin). In the isotropic case S'+_ couples to a subset of the 
states probed by Szz while in the anisotropic case the two functions probe different sets 
of states. 

In Fig. 2a we show the spin excitations that contribute to SzziQ, i^) for = 16 and 
for A = 1. As mentioned above the spectral weight of the spin excitations decreases 
as we move away from the lower threshold. The states contained in the region defined 
by e/(g) and eu{q) have a considerable weight while those at higher energies have a 
much smaller weight 0. The spectral weight of these higher states will vanish in the 
thermodynamic limit. Also other states contained in the continuum have very small 
weights. The continuum is therefore well defined by the set of states with highest 
spectral weight. In Fig. 2b we show the structure function for the Heisenberg chain. 
The delta functions at the excitation energies have been given a finite width both in 
frequency and momentum for better visualization. 

In Fig. 3 we show the states with non-vanishing spectral weight for the longitudinal 
and the transverse structure functions for the XY chain. The spectral weight is uniform 
inside the continuum defined by eqs. (6). The spectrum of the transverse excitations 
is however different. The lower spinon dispersion is well described by the single-spinon 
dispersion; particularly close to g = vr the gap is already rather small for such a small 
system. At higher energies the spectral weight is considerably spread. 

As we introduce the next-nearest-neighbor interaction the spectrum remains gapless 
for all A if j is small. In the intermediate coupling regime (j ~ 1) the system becomes 
gapped. In Fig. 4 we show the a) spectral weight and b) the structure function for 
j = 1 in the isotropic case. The states with 5*2 = and momenta = 0, vr are nearly 
degenerate |3^. The next excited state is a spin triplet with momentum k = 7t/2 
which in the thermodynamic limit will have a gap to the groundstate. Accordingly the 
spectral weight shows a gap with a continuum above it indicative of massive spinous. 
In the anisotropic case the two lowest states are the same as in the isotropic case but 
the next excited state is a Sz = 1, k = 7t/2 state The next state is a 5*2 = 0, /c = vr 
singlet. The low value of the gap signals the near level crossing that for = 16 occurs 
around j = 1.2 leading to a gapless regime ||3^. In Fig. 5 we show the spectral weights 
and the structure functions for the longitudinal and the transverse spin excitations. In 
the case of Szz the lowest gap is at = vr, while for S+^ the lowest gap is at = n/2. It 
is also clear that the spectrum is quite sharp at = 7r/2 in the transverse spin function. 

As we increase j further the sharpness of the gapless transverse mode at k = 7r/2 
becomes stronger. In Fig. 6 we show the structure function for the isotropic case and 
in Fig. 7 the same function for the anisotropic case at j = 2. The Goldstone mode 
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predicted to occur in the anisotropic case for the transverse spin excitations clearly 
singles out. 

We also consider the finite-size dependence of the low-energy excitations for the 
longitudinal and the transverse spectral functions as a function of the system size using 
results from exact diagonalizations and the modified Lanczos method. We consider 
system sizes up to = 24. The results extrapolate to the single-spinon dispersion 
curve for the various values of A and j. In particular we consider the anisotropic case. 

In Fig. 8 we show the lowest energy states for a) = and b) 5*^ = 1 as a 
function of momentum for the XY case for j = 2. The spectrum clearly shows the 
doubling of the lattice cell with a significant low energy mode at g = 7r/2, particularly 
in the transverse correlations {S^ = 1) where once again as the system size increases the 
gapless nature of the spectrum is evident as j grows (weak interchain couplings) . 

In the single-chain case the longitudinal spectrum can be obtained considering 
the two-spinon curves (assuming non-interacting spinous) via the usual procedure 
E{q) = e{ki) + e(/c2) where q = ki + k2, E{q) is the two-spinon curve and e(fc) is 
the single-spinon dispersion curve. These two limiting curves define the region of the 
continuum spectrum. The transverse excitations in the anisotropic case probe however 
a different set of states as can be seen for example from the single-chain XY case shown 
in Fig. 3c. The high energy part of the spectrum shows that the interactions between 
the spinous can not be ignored (remember that in the Jordan- Wigner transformation 
from a XY chain to spinless fermions, the transverse correlations involve the presence 
of strings). Therefore the two-spinon rule requires a proper treatment of the spinon 
interactions. The zig-zag ladder case is still more involved particularly for the transverse 
excitations. In any case the continuum is clearly visible. 



5. Conclusions 



The finite-energy and finite-momentum structure function provides a direct way of 
analysing the excitation spectrum. Previously the structure factor was analysed in 
the isotropic case ||26|, ^ ^ . The structure factor is obtained integrating the structure 
function over frequency at a fixed momentum. In the isotropic case the peak in the 
structure factor shifts from g = 7rtog = 7r/2 when the spiral phase is reached as the 
Majumdar- Ghosh point is crossed. In the anisotropic case on the other hand it has been 
predicted that in the limit of very weak interchain coupling an incommensurate gapless 
chiral phase should be observed [B^. However, for finite systems it is difficult to detect 



the incommensurability since the shift from commensurability is very small pO |. 

In this work we have analysed the structure function itself for the zig-zag ladder as 
a function of anisotropy (A) and the interchain coupling (j). The results show that in 
general the excitations are gapless or gapped spinous that have to be created in pairs as 
for the single-chain case. In the XY case (A = 0) as j grows it is clear that a gapless 
mode in the transverse excitations arises in agreement with previous results from field 
theory |2^ and with previous numerical results ||3^ obtained analysing the stiffness and 
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the level crossings. In the isotropic case (A = 1) the spectrum is a continuum of massive 
spinous. These results may be relevant to understand the ladder limit in the context of 
the recent experimental results on the two-dimensional system Cs2CuCU ||l|- 

After this work was completed we became aware of a preprint where using 
Mobius boundary conditions it is shown that in the isotropic case in the strong and 
intermediate coupling cases the spectrum may be described by a continuum resulting 
from two-spinon scattering, in agreement with the general conclusions of our paper. 
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Figure Captions 



Fig. 1- Excitation energies from the exact diagonalization oi a, N — 16 chain as a 
function of momentum for Sz — and Sz — 1 for the Heisenberg chain and the XY 
chain. 

Fig. 2- a) Excitation energies of the states that contribute to the longitudinal spectral 
weight for the Heisenberg chain. The solid lines arc the exact Bethe ansatz results for 
the thermodynamic limit. The numerical results arc obtained via exact diagonalization 
of a iV = 16 system. The color code of the points is the following: the points correspond 
to states with a spectral weight that is i) M{q,uj) > 1 (black), ii) 1 > M{q,u) > 0.1 
(red) and iii) 0.1 > M{q,uj) > 0.01 (green). The same color codes are used in the 
remaining figures. In Fig. 2b we show the structure function. The momentum is shown 
in units of tt. The vertical scale is in arbitrary units. 

Fig. 3- Excitation energies of the states that contribute to the a) longitudinal and c) 
transverse spectral weight of the XY chain. The solid lines in Fig. 3a are the exact 
Bethe ansatz results. We also show the structure function for the b) longitudinal and 
d) transverse excitations. 

Fig. 4- a) Excitation energies obtained from the exact diagonahzation of a A?" = 16 
ladder in the isotropic case and j — 1 that contribute to the spectral weight and b) the 
structure function. 

Fig. 5- a) Excitation energies of the states that contribute to the a) longitudinal and c) 
transverse spectral weight in the anisotropic case and j = 1. We also show the structure 
function for the b) longitudinal and d) transverse excitations. 

Fig. 6- Structure function of a = 16 ladder in the isotropic case and j — 2. 

Fig. 7- Structure function for the a) longitudinal and b) transverse excitations in the 
anisotropic case and j — 2. 

Fig. 8- Lowest energy branch for A = and j = 2 for a) 5*^ = and h) Sz = 1 for 
different system sizes. 
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